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ABSTRACT 

We present UV luminosity functions (LFs) at 1500 A derived from the HST Early Release Science 
WFC3/UVIS data acquired over - 50 arcmin^ of the GOODS-Soutli field. The LFs are determined 
over the entire redshift range z = 0.75 — 2.5 using two methods, similar to those used at higher 
redshifts for Lyman Break Galaxies (LBGs): (1) 13-band UV+optical+NIR photometric redshifts to 
study galaxies in the range z = 0.5 — 2 in three bins of dz = 0.5, and (2) dropout samples in three 
redshift windows centered atz^l.5,z^l.9, and z ^ 2.5. The characteristic luminosity dims by 1.5 
mag from z = 2.5 to z = 0.75, consistent with earlier work. However, the other Schechter function 
parameters, the faint-end slope and the number density, arc found to be remarkably constant over the 
range z = 0.75 — 2.5. Using these LF determinations we find the UV luminosity density to increase by 
~ 1.4 dex according to (1 + z)'^-^^^'^-^^ from z '-^ to its peak at z ^ 2.5. Strikingly, the inferred faint- 
end slopes for our LFs are all steeper than a — —1.5, in agreement with higher-redshift LBG studies. 
Since the faint-end slope in the local universe is found to be much flatter with a ~ —1-2, this poses 
the question as to when and how the expected flattening occurs. Despite relatively large uncertainties, 
our data suggest a ~ —1.7 at least down to z ~ 1. These new results from such a shallow early dataset 
demonstrate very clearly the remarkable potential of WFC3/UVIS for the thorough characterization 
of galaxy evolution over the full redshift range z^0.5toz'^3. 

Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: luminosity function, mass 
function 



1. INTRODUCTION 

The UV Luminosity Function (LF) of galaxies is a key 
diagnostic for establishing the contributions of galaxies 
of different luminosities and masses to the cosmic star- 
formation rate (SFR) density. In addition, in combina- 
tion with other measures, like the galaxy mass function, 
it also provides essential clues as to the causes for the 
slow down in star-formation of galaxies in more recent 
times since z '^ 2. 

The UV LF has been studied very extensively at red- 
shifts 2 < z < 6, but at lower redshifts (z < 2), 
similar deep, high resolution data has not been avail- 
able. To date the z < 2 UV LF has been measured 
from GALEX data using the dropout technique (e.g. 
Burgarella et al. 2006; Ly et al. 2009; Haberzettl et al. 
2009), and also using spectroscopic or photometric 
redshifts (e.g. Arnouts et al. 2005; Wydcr et al. 2005; 
Budavari et al. 2005). The main limitations of the 
GALEX results are (1) the very wide point-spread func- 
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tion (PSF; ~ 5" FWHM), which may cause blending of 
the UV light from several sources and makes the iden- 
tification of optical counterparts relatively difficult, and 
(2) the limited depth (^ 25 mag). 

With the installation of the Wide-Field Camera 3 
(WFC3) on the Hubble Space Telescope (HST), the ef- 
ficient coverage of the electromagnetic spectrum at high 
spatial resolution (<0'.'15) has been extended to the UV 
and to the IR. The new IR capability allowed us (and 
other groups) to perform the first robust measurements 
of the UV LF at z ~ 7 — 8 in the first-epoch data of the 
HUDF09 (e.g. Ocsch ct al. 2010b; Bouwens et al. 2010; 
McLure ct al. 2010; Bunker ct al. 2010). 

The high sensitivity UV imaging capability provided 
by the WFC3/UVIS channel similarly allows us to de- 
rive the UV LF at redshifts below z ^ 2.5, where the 
'--^1500 A UV flux from galaxies could not be probed ef- 
ficiently with previous high resolution cameras. While 
the depth of the current ERS WFC3/UVIS UV data 
(Windhorst et al. 2010) is less than that of the WFC3/IR 
data, the present depths still significantly exceed that 
from any previous UV data set. 

We use two different methods for our LF determina- 
tions: (1) 13-band UV-|-optical-t-NIR photometric red- 
shifts to study galaxies in the range z ^ 0.5 — 2, and 
(2) the UV dropout technique, applied to the UVIS filter 
set, which allows for the selection of star-forming galax- 
ies at z ^ 1.3 to z ~ 2.8. The latter dropout samples are 
similar to the recent selections of Hathi et al. (2010). 

In Section 2 we describe the data, the photometric red- 
shift estimates and the dropout selection. This is fol- 
lowed in Section 3 by our constraints on the UV LF be- 
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tween z = 0.5—2.8, and wc end with a short discussion on 
the evolution of the UV LF parameters and luminosity 
density in Section 4. 

Throughout this letter, we adopt $7^/ = 0.3, 51a = 
0.7, iJo = 70 kms^^Mpc-^ i.e. h = 0.7. Magnitudes 
are given in the AB system (Oke & Gunn 1983). 

2. DATA AND SAMPLE SELECTION 
2.1. ERS Data Set 

The Early Release Science (ERS) UV data cover 
a total of ^50 arcmin^ of the GOODS-South fields 
(Giavalisco et al. 2004) in 4 x 2 pointings of WFC3/UVIS 
in three filters {UV225 - F225W, UV275 ^ F275W, Upe 
- F336W). Each pointing was covered with five orbits, 
equivalent to 2778 s, 5688 s, and 5688 s integrations 
in C/V225, UV275, and f/336, respectively. The observa- 
tions in each band were split into three (t/336) or four 
(t/V225/275) exposures to facilitate cosmic-ray (and bad 
pixel) removal (see Windhorst et al. 2010, for more in- 
formation on the data). The RMS maps were scaled to 
the proper noise level, as measured in circular apertures 
appropriate for the galaxies under study. Using circular 
apertures of 0'.'2 radius, the measured 5cr point source 
sensitivities in all three filters are 26.8-26.9 mag. 

After masking areas with bad sampling and large resid- 
uals from cosmic ray hits the final area used in this 
study with both ACS and WFC3/UVIS coverage is ~ 47 
arcmin^ . 

The reduced UVIS images were registered to the 
GOODS ACS B-band frames (version 2; M. Giavahsco 
and the GOODS Team, 2010, in preparation^) and driz- 
zled to the same pixel scale of 0'.'03. The most recent 
in-flight geometric distortion solutions (02/2010) were 
essential for obtaining good alignment. Sources were 
detected in the ACS B-band image using SExtractor 
(Berlin & Arnouts 1996), resulting in 6683 detections 
down to 28 mag. Colors were measured in dual image 
mode in isophotal apertures to maximize the S/N. 

Lastly, our catalog was matched against the pub- 
lic GOODS catalogs and candidate stars flagged. At 
Z850 < 25.5 this identification is based upon their size 
(half-light radii < 0'.'09) and the SExtractor stellarity 
parameter (> 0.8). Faintward of 25.5, the latter param- 
eter is no longer reliable, and therefore flagging sources 
with half-light radii < 0'.'09 was done if the V — i vs. 
i — z colors of sources were consistent with the stellar 
sequence (within 0.15 mag). From a visual inspection of 
the sample of galaxies entering in the LF estimates this 
has proven to be a very efficient and complete procedure 
to identify stars. 

2.2. Photometric Redshift Selections 

In order to derive accurate photometric redshifts, we 
supplement the high resolution ACS and WFC3/UVIS 
data with ground-based IR and with Spitzer data, using 
the GOODS-MUSIC catalog (v2) of Santini et al. (2009). 
In particular, the GOODS-MUSIC photometry used here 
is VIMOS [/, ACS 5435/^606/^775/^850, ISAAC J,ff, A", 
and Spitzer IRAC 3.6 and 4.5. The GOODS-MUSIC 
catalog is 90% complete down to 2350 < 26 mag. Match- 
ing this catalog to our B-selected sources, we find 4165 

^ http://archive.stsci.edu/pub/hlsp/goods/v2/ 



matches (using a 0'.'2 error tolerance). The UVIS isopho- 
tal fluxes are scaled to total fluxes using the ratio of 
B-Band ISO flux to B-Band AUTO flux. 

Thus, our basic photometry catalog for photo- 2: esti- 
mates consists of 7 filter HST data, complemented with 
4 PSF-matched ground based filters, and the two shortest 
wavelength IRAC bands. The IRAC 5.8 and 8.0 bands 
were omitted, since they are dominated by dust and PAH 
emission for lower redshift galaxies. Even if galaxies are 
generally red in B — z, we limit our analysis to i?435 < 26 
mag to avoid incompleteness corrections due to the ^850- 
band limit. We did not use the WFC3/IR photometry, 
since it is not available for all sources to keep our analysis 
as uniform as possible (and our tests show that our se- 
lections do not change significantly when these data are 
included) . 

Photometric redshifts were estimated from spec- 
tral energy distribution (SED) fitting with ZEBRA 
(Feldmann et al. 2006). The adopted SED set is based 
on Bruzual & Chariot (2003) composite stellar popula- 
tions. In particular, we computed 15 median templates 
as a function of rest-frame u — J color for galaxies in the 
zCOSMOS survey (Lilly et al. 2007; Lilly et al. 2009). 
Subsequently, these templates were corrected for system- 
atic offsets with respect to the COSMOS photometry 
(Capak ct al. 2007) using ZEBRA, and the major emis- 
sion lines (Lya, Ha, H/3, OH, and OIII) were added using 
the Kennicutt (1998) relations between UV luminosity 
and emission line strengths. Finally, dust was added to 
this set of templates using the Calzetti et al. (2000) dust 
law with E{B — V) = Q — 0.5, resulting in a final set of 
2052 templates. 

We used this template set to derive photometric red- 
shifts for our combined catalog of UVIS and GOODS- 
MUSIC photometry. We compared our estimates with 
spectroscopic redshift measurements compiled in the 
GOODS-MUSIC catalog. From a total of 462 sources 
with reliable spectroscopy (363 at z < 1.5), we find an 
accuracy of (Tz = 0.037 x (1-l-z), with only 6.4% of sources 
having redshift errors larger than |A(z)| > 0.15 x (1-l-z). 
Our photometric redshift estimates are therefore suffi- 
ciently precise to derive LFs in bins of dz = 0.5. 

2.3. UV Dropout Selection 

Most of the UV LFs at z > 3 are based on an LBG 
selection (e.g. Stcidcl et al. 1995). The same technique 
can also be applied to the WFC3/UVIS filters, allowing 
the selection of star-forming galaxies as UV dropouts. At 
z < 3, the Lyman break seen in galaxies is primarily the 
result of neutral hydrogen within the sources themselves, 
in significant contrast to the situation at z > 3, where 
the break is also due to Lyman-series absorption from 
the IGM. 

Our selection criteria are based on standard color-color 
diagrams, which are shown in Figure 1. In particular, 
the selection regions are chosen to identify an unbiased 
sample of star-forming galaxies with dust extinction less 
than E{B — V) ^ 0.3 — 0.4 guarding against lower redshift 
interlopers. 

Specifically, the adopted selection criteria for UV225 
dropouts are: 
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Fig. 1. — The color selection of UV-dropouts at redshifts z ~ 1.5, 1.9, and 2.5 (from left to right). Blue lines represent tracks of star- 
forming galaxies with dust obscuration of E{li — V) = 0, 0.15, 0.3. The green dashed lines are the expected colors of lower redshift galaxies 
(Coleman et al. 1980). These tracks extend to 2 < 1.1, z < 1.5, and z <1 for the different figures respectively. The gray shaded regions 
are the adopted selection regions for dropout galaxies. The orange points are stars from the Pickles (1998) library. Stars are removed from 
our samples using a combination of the stellarity parameter and {V — i) vs. (i — z) colors, as described in the text. The light gray dots 
represent all galaxies in the catalog down to the respective magnitude limits adopted in our LF determinations (see section 3.2), and black 
squares indicate our candidates. Fluxes with a detection significance below Icr have been replaced with Icr upper limits. Typical errorbars 
are shown in the upper right. 



-0.5 < (C/F275 - C/336) < 1.4 A 

S/N{UV275) > 5 

This selects galaxies between z ~ 1.3 — 1.7, and results 
in 72 candidates down to C/336 = 26.5. 

Similarly, for t/V275-dropouts, the selection criteria 
are: 

UV275 - U336 > 1 A 
UV275 - f/336 > 2.2(C/336 - S435) - 0.42 A 

-0.5 < (C/336 -S435)< 1.1 A 

S/NiU336)>5 A 5/7V(t/F225) < 2 

selecting galaxies between z ~ 1.7 — 2.1, with 131 candi- 
dates down to i?435 = 26.5. 
And finally, C/336-dropouts are selected with: 

C^336 - B435 > 1 A 
C^336 — ^435 > 2(5435 - ^606 ) A 

-0.5 < (B435 - V^eoe) < 1 A 

5/7V(S435)>5 A S/N{UV225/275)<2 

returning 434 sources down to Vgoe = 26.5 with expected 
redshifts in the range z ~ 2.2 — 2.8. When applying these 
criteria, all fluxes below their Icr limits (in the dropout 
band) are replaced with the respective Icr upper limits. 

For the C/336-dropouts, the stellar sequence runs ex- 
actly through the selection window (see Figure 1). How- 
ever, stars are removed prior to candidate selection using 
the same identification scheme as for the photo-z sam- 
ple, resulting in a total of 40 excluded stars. From a 
final visual inspection of all sources no additional stars 
were found. A small number of clear diffraction spikes 
and spurious detections on the wings of large low-redshift 
galaxies were removed by hand. 

Using our photometric redshifts, we estimate the low- 
redshift interloper fraction of these UV-dropout samples 
to be very small. In particular, only 7% of the sources 
in the t/336 dropout sample have Zphot < 1.75. Similarly, 
only 4% and 8% of the sources in our C/V275- and UV225- 
dropout samples are interlopers, i.e., have Zphot < 1-25 
and Zphot < 1, respectively. Only sources down to the 



limiting magnitudes of our UV LFs are considered in 
this estimate (see section 3.2). 

3. RESULTS 
3.1. Photometric Redshift LF Results 

The UV luminosity functions at 1500 A in the range 
z ^ 0.5 — 2 are measured in three bins of photometric 
redshift of width dz = 0.5. The absolute magnitudes are 
computed from the UV275 , C/336 and i?435 photometry for 
sources in the redshift range z = 0.5 — 1, z = 1 — 1.5, and 
z = 1.5 — 2, respectively. K-corrections are derived from 
the appropriate best-fit templates. 

The UV luminosity functions of galaxies are estimated 
using the standard Knax estimator (Schmidt 1968) such 
that 

where the sum runs over all galaxies in the given redshift 
and absolute magnitude bin. We only include galaxies 
with > 5cr detections in the corresponding filter, which 
results in the need for a completeness correction factor C. 
Following Oesch et al. (2007, 2010b), this is determined 
from Monte-Carlo simulations in which we add artificial 
galaxies to the images and detect them with the same 
SExtractor parameters as used for the original B-band 
catalog. The completeness is > 97% to < 25 mag in the 
UVIS filters, after which it drops to 50% at 25.8 mag, 
which is our faint-end limit. 

Figure 2 shows the resulting luminosity functions. The 
best-fit Schechter parameters (Schechter 1976) are tabu- 
lated in Table 1. While we sample the faint end rather 
well, with data that reaches ^ 2 mag fainter than M^,, 
the relatively small area covered by the ERS data results 
in some uncertainty on the cutoff at bright magnitudes. 
Nevertheless, our results are in good agreement with the 
LFs determined from GALEX by Arnouts et al. (2005) 
and effectively validate those results in data with ^ 50 
times higher spatial resolution. 
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Fig. 2. — Left — The 1500 A LF in the three different bins of photometric redshift. The blue data points with errorbars are the step-wise 
LFs determined with the Vinax method, and the soUd blue Unes represent our best-fit Schechter LFs to the step-wise LFs. The best-fit 
parameters are presented in Table 1. The black dashed line is the UV LF in the local universe from Wyder et al. (2005), while the gray 
dashed line is the UV LF between z ~ 1.9 - 2.7 from Roddy & Steidel (2009). Right— The 1500 A LF derived from the LBG samples 
together with their best-fit Schechter LFs. From top to bottom: UV225-, UV275-, and f733g-dropout results. Again, the local z ~ and the 
z ~ 2.3 LFs arc shown for reference as black and gray dashed lines respectively. The insets show the expected redshift distributions of the 
dropout galaxies. This is estimated from the effective volumes and the best-fit Schechter LFs. 



TABLE 1 

UV Luminosity Function Parameters^ 



z 


No. 


of Sources 


logio'/'* [Mpc ^'mag i] 


M, 


a 




Photometric redshift samples 


0.5-1.0 
1.0-1.5 
1.5-2.0 




284 
284 
288 




-2.52 ±0.31 
-2.90 ±0.25 
-2.63 ±0.23 


-19.17 ±0.51 
-20.08 ± 0.36 
-20.17 ±0.34 


-1.52 ±0.25 
-1.84 ±0.15 
-1.60 ±0.21 


UV-dropout samples 


1.5 (Az = 0.4) 
1.9 (Az = 0.4) 
2.5 (Az = 0.6) 




60 
99 
403 




-2.64 ±0.33 
-2.66 ±0.36 
-2.49 ±0.12 


-19.82 ±0.51 
-20.16 ±0.52 
-20.69 ±0.17 


-1.46 ±0.54 
-1.60 ±0.51 
-1.73 ±0.11 


t Derived from t. 


'2 fits to the 6 


tep 


-wise LFs. 







3.2. Dropout LF Results 

The UV LF of the UVIS dropout sample is computed 
using standard techniques adopted in LBG studies. In 
particular, the step-wise LFs are derived using 

(/.(MK, z])dM = NjV,ffim,) 

where N^ are the number of galaxies observed in the mag- 
nitude bin i, and the effective volume is given by: 

f°° dV 

Vesinii) ^ S{mi,z)C{mi)—dz 



The selection efficiency S and completeness C are again 
estimated from simulations of artificial galaxies inserted 
into the observed images. The inserted galaxies are 
chosen to have a log-normal size distribution, using 
an extrapolation of the higher redshift LBG size scal- 
ing of (1 -I- z)^^ (Oesch et al. 2010a; Bouwens et al. 
2004; Ferguson et al. 2004). The colors are set accord- 
ing to the z '^ 2.5 UV continuum slopes derived in 
Bouwens et al. (2009) (see also Reddy et al. 2008), with 
additional IGM absorption according to Madau (1995). 
Due to the lower sensitivity of the UV data with re- 
spect to the optical, red colors across the continuum 
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breaks cannot be accurately measured for galaxies fainter 
than B ^ 26. Therefore only sources with V < 26.3, 
B < 25.5, and C/336 < 25.4 are used to compute the 
LFs, ensuring that the selection efficiency is still larger 
than 50% for C/336-, UV275-, and C/V225-dropouts, re- 
spectively. This results in 403, 99, and 60 sources in the 
three different dropout samples. 

The UV LFs for the UVIS dropouts are shown in the 
right panels of Figure 2. Again, a significant increase at 
the bright end of the UV LF is seen from lower to higher 
redshift. Due to the modest size of the present UV225- 
and C/V275-dropout samples, the Schechter parameters 
we derive at z ^ 1.5 and z ^ 1.9 are still somewhat 
uncertain (see Table 1). The z '^ 2.5 [/336-dropout LF, 
by contrast, is much better constrained from our data, 
extending '^ 2 mag fainter than Af* and with a faint- 
end slope a established to within ±0.11. Our z '^ 2.5 
LF is in excellent agreement with the z ^ 2.3 UV LF of 
Reddy & Steidel (2009). 

Note that due to the limited area of the present survey, 
our results are dominated by cosmic variance, which we 
estimate to be about 16 - 20% (Trenti & StiavelH 2008). 

4. DISCUSSION 

The evolution of the Schechter parameters of recent 
UV LFs in the redshift range 2 ~ — 8 is shown in 
Figure 3. Our results cover the whole 4.3 Gyr period 
from z ^ 0.75 to z ^ 2.5. Interestingly, in the region of 
overlap, the LFs we derive from our photometric redshifts 
and UV dropout selections are in very good agreement. 
This suggests that our UV dropout samples are well- 
defined and reasonably complete. 

In Figure 3, our results are compared to several 
Schechter function estimates from the literature. Within 
the uncertainties, 0» is found to be fairly constant over 
the entire redshift range, although there is a weak trend 
towards higher (p^ at lower redshifts. 

The most noteworthy trend is the monotonic fading of 
the characteristic luminosity by a factor of ~ 16 from 
z '^ 3 (Af* fades from —21 at z ^ 3 to —18 at z ^ 0). 
At z > 4, M* turns over and decreases monotonically 
again towards higher redshift, such that M^,{z = 8) is 
essentially equal to Af*(z ^ 1). 

It is very interesting that all our faint-end slope mea- 
surements are steeper than a ^ —1.5, even at z ^ 
0.5 — 1, in agreement with the earlier measurements of 
Arnouts et al. (2005). Despite relatively large crrorbars, 
the faint-end slope a is clearly seen to transition from 
a < —1.5 at z > 0.5 — 2 to a flatter a — —1.2 in the local 
universe (but see Treycr et al. 1998, who find a steeper 
slope locally). With the present data we cannot deter- 
mine whether this transition occurs relatively abruptly 
at low redshift, or smoothly from z ~ 2 to z '^ 0. 

Finally, in Figure 4 we show the cosmic 1500 A lumi- 
nosity density pL (uncorrected for dust extinction) esti- 
mated from the Schechter function parameters. The lu- 
minosity density is shown integrated to L > 0.1L*^3 and 
also integrated to zero luminosity. Clearly, the evolution 
of PL is mainly driven by Af*. From the peak at z '^ 2.5 
the total UV luminosity drops by 1.5 dex to z '^ 0. This 
decrease is well described hy pl{z) oc (1 -|-z)"^ , for which 
we find /3l(z < 2.5) = 2.58 ± 0.15, in good agreement 
with earlier measurements (e.g. Schiminovich et al. 2005; 
Tresse et al. 2007). The evolution at higher redshift fol- 
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Fig. 3. — Top Three Panels — The redshift evolution of UV 
(A = 1500 — 1700 A) LF parameters <^«, M*, and a over 
the redshift range z ~ — 8. The present determinations 
are shown with large filled circles. Of particular note, the un- 
certainties on the Schechter function parameters are frequently 
much smaller at 2: > 3 than at lower redshifts. This is a 
direct consequence of having very deep HST optical surveys. 
WFC3/UVIS should enable similarly-deep future observations in 
the UV at z < 3. Also shown arc some notable determinations 
from the literature (Arnouts et al. 2005; Reddy & Steidel 2009; 
Bouwens et al. 2007, 2008; Bouwens et al. 2010; Oesch et al. 2007, 
2010b; Sawicki & Thompson 2006; Yoshida et al. 2006; Iwata et al. 
2007; McLure et al. 2009; Ouchi et al. 2009). Despite some scatter, 
most faint-end slope determinations at 2 > 0.5 are very steep, i.e. 
a ~ —1.7. Bottom — The 68% error contours on M* and a over 
the redshift range z ~ 0.75 — 2.5. The contour at 2 ~ 2.5 is from 
our C/336-dropout sample, the remaining ones are derived from the 
photometric redshifts. The lower redshift dropout samples result 
in considerably larger errors and are omitted for clarity. 

lows (3l{z > 2.5) = -1.82±0.37. z ~ 2.5 thus marks the 
high point in the production of the cosmic UV luminos- 
ity and is clearly an important epoch for exploring the 
properties of galaxies (e.g. see Cameron et al. 2010, for 
a new technique for selecting galaxies in this range with 
simple color criteria). 
The remarkable potential of the HST WFC3/UVIS 
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camera is clear from these new results. This ERS dataset 
with its initial rather shallow observations has demon- 
strated already that WFC3/UVIS is a powerful tool for 
exploring galaxy evolution over the redshift range z '^ 0.5 
to z ~ 3. Deeper data will provide key insights into the 
decline of the cosmic SFRD since z ~ 2, particularly at 
the lowest luminosities. For example, a steep LF with 
a ^ —1.7 suggests that star- formation efficiencies and 
feedback processes are only weakly dependent on halo 
mass (e.g. Trcnti et al. 2010). 



Fig. 4.— The redshift evolution of the UV (A ~ 1500 A) lumi- 
nosity density. In the upper panel, the luminosity functions are 
integrated down to a fixed luminosity of 0.1L*_3 (L* as derived 
at z = 3 by Steidel et al. 1999), while the lower panel shows the 
total luminosity density. The dark and light blue points are de- 
rived from our photo-2 and the dropout samples respectively. The 
gray points are based on the literature values shown in Figure 3. 
The dashed curves are fits to the data at z < 2.5 and z > 2.5 of 
the form pj^ oc (1 -|- z) "'^ . Errorbars are computed from Monte- 
Carlo simulations, not accounting for the covariance between the 
Schechter parameters. They are thus expected to overestimate the 
real errors. Note, however, that a was constrained to a > —2 in 
order to keep p^°' finite. This only significantly affected 5 data 
points, which arc marked with black circles. 
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